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We consider a nonequilibrium reaction-diffusion model on a finite one dimensional lattice with 
bulk and boundary dynamics inspired by Glauber dynamics of the Ising model. We show that the 
model has a rich algebraic structure that we use to calculate its properties. In particular, we show 
that the Markov dynamics for a system of a given size can be embedded in the dynamics of systems 
of higher sizes. This remark leads us to devise a technique we call the transfer matrix Ansatz 
that allows us to determine the steady state distribution and correlation functions. Furthermore, 
we show that the disorder variables satisfy very simple properties and we give a conjecture for 
the characteristic polynomial of Markov matrices. Lastly, we compare the transfer matrix Ansatz 
used here with the matrix product representation of the steady state of one- dimensional stochastic 
models. 
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I. INTRODUCTION 



The study of systems far from equilibrium has been greatly helped by the discovery of exactly solvable models, 
because explicit computations for these models provide us with excellent testing grounds for general hypotheses about 
nonequilibrium statistical mechanics [J, Q . 

One important difference between equilibrium and nonequilibrium behaviour is encoded in the detailed balance con- 
dition. This condition states that at equilibrium the total transition rate between two arbitrary micro-configurations 
vanishes identically Conversely, nonequilibrium steady states usually break detailed balance, which results in 
the existence of current loops in the configuration space of the system, leading to a non-zero macroscopic (physical) 
current that transports matter, momentum or energy from one region of the system to another. Because of this 
current, the boundaries of the system can affect its bulk and the modification of the boundary conditions through a 
control parameter can induce dynamical phase transitions even in one dimensional systems Such sensitiveness to 
the boundary conditions is well demonstrated by exact results obtained for the asymmetric exclusion process (ASEP) 
in one dimension, which is one of the simplest examples of a driven lattice gas and one of the most exhaustively 
investigated interacting particle systems d, 0, 0| ■ For the ASEP on a periodic ring the steady state is uniform and 
all configurations have the same stationary weight; in contrast, for the ASEP on a finite lattice with open boundary 
conditions (that allow injection and removal of particles at the end sites) the measure is non-uniform and in the limit 
of large sizes, the system can exist in three different phases: maximal current, low density and high density (the last 
two phases being separated by a line of shocks). The exact expression of the stationary measure valid for any system 
size was first derived in Q, introducing a method now called the matrix product representation that has become an 
important technique for one dimensional interacting particle processes. In particular, this matrix representation has 
an interpretation in terms of discrete lattice paths that leads to exact combinatorial results for finite size systems (see 
for an exhaustive and recent review). Thanks to this matrix product method and to more standard integrability 
techniques such as the Bethe Ansatz , an plethora of results have been derived for the ASEP [ll[ . 

Here, we study a reaction-diffusion model on a finite lattice of L sites with open boundaries in which hard-core 
particles perform asymmetric jumps and can undergo pairwise annihilation. This model is inspired by the Glauber 
dynamics [T^ for the Ising model, in which the elementary excitations (or particles) are not the individual spins but 
rather the domain walls between sets of opposite spins. More precisely, we shall consider the totally asymmetric 
Glauber dynamics where each spin changes its orientation with a certain probability based strictly on the spin to its 
left. Hence, the domain walls, which are represented by particles, move only to the right and, if two of them collide, 
they annihilate each other. Reaction-diffusion models have been thoroughly studied in the nonequilibrium statistical 
physics literature (see e.g. [l^ for similar studies on the real line and [H, [l^, [13, [H, [3 '^^ the lattice) and 
their relation to non-Hermitian spin chains has been established [l^, |2l[. In particular, it has been shown that for 
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certain values of the reaction-rates, they reduce to free fermion models that can be solved by using Jordan- Wigner 
fermionization techniques [13, [H, 0, [U H^, [13, [13 • This method was applied to periodic boundary conditions where 
translation invariance allows the use of Fourier transform to diagonalize the free fermion Hamiltonian. More recently, 
similar techniques were used by Farago ^2§\ , and Farago and Pitard 30] to calculate large deviation functions of the 
time integrated injected power when a single spin is allowed to perform Poissonian flips. 

We shall study the effect of open boundaries conditions on the long time dynamics of the model. At the boundary 
sites of the lattice, we allow creation and destruction of particles in such a manner that the boundary dynamics is 
compatible with the bulk dynamics. Here, the Jordan- Wigner transformation introduces non-local terms in the free 
fermion Hamiltonian and the lack of translation invariance precludes the use of the Fourier transform. 

We shall define a one-parameter family of models with a nontrivial steady state that generalizes the totally asym- 
metric exclusion process (TASEP) on a finite lattice, first solved in Q. For a special value of the parameter (which 
corresponds to the free fermion point), we calculate the exact noncquilibrium steady state measure by using a recur- 
sion between systems of different sizes. Here, the key to the solution is not a matrix product representation but a 
different approach: using a linear transform between configuration spaces of different dimensions, we shall prove that 
the Markov matrix of the system of size L can be embedded in the Markov matrix of the system of size L + 1. The 
existence of such a rectangular (thus, non-invertible) similarity matrix, that will be called a 'transfer matrix Ansatz', 
will allow us to derive exact combinatorial expressions for the local density and for correlation functions of the model, 
and to extract their asymptotic behaviour in the infinite system size limit. 

The organization of the paper is as follows: In Section |TT1 we introduce the model. In Section Hill we describe the 
transfer matrix Ansatz in general terms and apply it to our system. This leads us to a closed form expression for the 
'partition function' that normalizes the stationary probabilities and we derive a formula for the joint occupation of 
the first k sites for all k. In Section ITVl we define the disorder variables and calculate the associated one-point and 
two-point functions. We use these variables to calculate the density and the rate of evaporation; we also prove some 
general properties of higher correlation functions. In Section |Vl we conjecture some properties of the spectrum of the 
Markov matrices. Section IVII is devoted to concluding remarks and some open problems. In Appendix [21 we write 
explicit recursion relations for the steady state probabilities that allow us to compare the transfer matrix Ansatz with 
the matrix product representation. 



II. THE MODEL 



We consider a nonequilibrium system on a finite lattice with L sites labelled from 1 to L. We denote the boundary 
of the domain wall, that is the boundary between oppositely charged spins, by a particle using the standard notation 
1. Empty sites are denoted by 0. The evolution rule in the bulk, biased Glauber dynamics, is thus given by 

10 01 with rate 1 , 

1) 

11 -> 00 with rate A. ^ ^ 



We remark that the first rule represents the movement of the domain wall to the right (eg. + +| ^ + + +|^) 

and the second, the annihilation of two domain walls (eg. + + | ^ |h > + + ++)• 

The evolution of the first site is given by 

0^1 with rate a, , , 

(2) 

1^0 with rate aA. 



Particles can exit from the last site according to 

1-^0 with rate /3. (3) 

The rules ([2]) and ([3]) were constructed by considering the finite lattice to be a 'window' of the infinite one-dimensional 
lattice. Suppose that there is a virtual site labelled to the left of the first site and a virtual site labelled L -I- 1 to 
the right of the last site. The left boundary conditions are deduced from the bulk rules (jll by looking at the second 
component of the bond 0—1 and the right boundary conditions are obtained by looking at the first component of the 
bond L — (L + 1); a and /3 are free control parameters. 

The rules ([1]), ^ and ^ fully define the noncquilibrium model by allowing one to construct explicitly its Markov 
matrix of dimension 2^ . Notice that we use terminology directly from the original TASEP problem for the bulk and 
boundary rates and that we have ensured the consistency of the bulk and the boundary rates by fixing the rates 
appropriately. Indeed, the value A = corresponds exactly to the TASEP solved in [1]. 
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The model can be discussed using an approximate mean field argument. The evolution equation for the density at 
site z > 1 is given by 

^(7?,) = ~ (A + 1)7^,;)) - (^.(1 - (A - , (4) 

where rji represents the occupation variable at site i. For A = 1, these equations become simpler because the state of 
a site depends only on the preceding sites. Then, using the mean-field assumption, we obtain the following recursion 
valid in the stationary state: 

im) - (5) 
The stationary density at the first site can be obtained exactly by writing the equation 

0=4('/i> 2r7,)>-(?7i), (6) 



dt 

and therefore one obtains the general formula 



1 + 2ka 



This implies that the density falls off like k^^ for large k. This mean- field result is wrong: in one-dimension, the 
actual exponent is —1/2. Note that the exit rate /3 did not enter the calculation. 

The fact that the value A = 1 plays a special role can readily be understood from the dynamical rules: indeed the 
exit rate from any site i is equal to 1 whether site z -I- 1 is occupied or not. Besides, A = 1 corresponds to the free 
fermion point of the associated spin chain (sec Section IVTl for a more detailed explanation of this fact). 

Remark 1 More generally, for A = 1, this model has the following important property. The correlation function 
{rji-^ . . . rjij) does not depend on the state of ?7i„+i simply because the exit rate from any configuration that contributes 
to this expectation value is equal to 1 whether site + 1 is occupied or not (by ([T|) ). In particular, for a system of 
size L, all correlation functions which depend on sites other than the last one are going to be strictly independent of 
[3 and moreover, any correlation function of the form {rji-^ . . . rji J) is independent of L as long as L > in + I. 

For the rest of the paper, we will take A = 1 as that is the only case for which we can derive explicit combinatorial 
formulae. 



III. SEMI-SIMILARITY BETWEEN MARKOV PROCESSES 



A. The Transfer Matrix Ansatz 



Let AIl denote the Markov matrix for a system with L sites. Typically the size of the matrix will be exponential 
in L. We first give a general definition of the "transfer matrix Ansatz" and then apply it to the specific case of the 
problem defined above. 

Let us consider a family of Markov processes defined by Markov matrices {M^} of increasing sizes (in most physical 
applications, L is the size of the system). We shall see that for certain systems there exists a natural embedding of 
Ml into Ml+i- 

Definition 1 We say that a family of Markov processes satisfies the Transfer Matrix Ansatz if there exist matrices 
Tl,l+i for all sizes L such that 

Ml+iTl,l+i = Tl,l+iMl . (8) 
We also impose that this equality is nontrivial in the sense that 



Ml+iTl.l+i ^ 0. 



(9) 
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The rectangular transfer matrices T^^^+i can be interpreted as semi- similarity transformations connecting Markov 
matrices of different sizes. Anotlier way to view the transfer matrices is that that the following diagram 



Mr 



Ml + 1 



L+1 



(10) 



commutes, where flL is the space of 2^ configurations of size L. 

We first explain the importance of the last condition. For the noncquilibrium systems that we are interested in, 
there is usually one unique steady state, which means the multiplicity of the zero eigenvalue is one. If \vl) is a non-zero 
vector in the kernel of Ml and = (1, 1, ... , 1), the matrix Vl,l+i = \vl+i){^l\ satisfies (O since the Markov 
matrices satisfy the condition {1l\Ml — 0. However, we have Ml+iVl.l+i = Vl,l+iMl = 0, and the condition Q 
is violated: therefore this trivial solution is excluded. 

The above definition leads immediately to some important observations. First we have 



- Tl.l+iMl\vl) = Ml+iTl.l+i\vl) 



(11) 



which, assuming Tl,l-i-i\vl) ^ 0, and taking into account the uniqueness of the steady state, allows us to define |fL+i) 
so that 



Tl,l+i\vl) = \vL+i)- 

The other important consequence is related to the eigenvalues oi M^. Let 
II. Then 



(12) 

be any eigenvector of Ml with eigenvalue 



IJ-Tl,l+i\w) ^Tl,l+iMl\w) ^ Ml+iTl,l+i\w) 



(13) 



which, again assuming Tl^l+i\w) ^ 0, means that Tl^l+i\w) is an eigenvector of Ml+i with the same eigenvalue 
fi. In other words, an eigenvector of Ml that is not in the kernel Tl^l+i is also an eigenvector of Ml+i with the 
same eigenvalue (if the matrix Ml is not diagonalizable, then we cannot immediately make a statement about the 
multiplicities). Conversely, if Ml and Ml+i have a common eigenvalue /it, then the rank-one rectangular matrix 
is a transfer matrix; but the image of the steady state \vl) of Ml by this matrix vanishes. Hence, 
\lJbL+i){lJ'L\ cannot be used to construct the steady state of Ml+i knowing \vl)- Therefore, in order to study the 
stationary state, we must look for transfer matrices that satisfy the additional condition: Tl^l+i\vl) 0. 



B. The case of the asymmetric annihilation process 

For the system we consider here, the Markov matrices are of size 2^. It is most convenient for us to take the 
naturally ordered basis of binary sequences of size L. For example, when L = 2, the ordered list is (00,01, 10, 11). 
The first important observation is that there is a recursion of order one among the Markov matrices. 



Theorem 1 Let a denote the matrix 

/n 1 \ 

(14) 



1 

1 



and \l denote the identity matrix of size 2 . Then 



AT ^ f Ml^i - a{a (g) 1l-2) oIl-i + {a (g> 1l-2) \ t.^ 

^ V alL-l ML-l-lL-l-a{<J<»lL-2) J ' ^ ' 

where Ml is written as a 2 x 2 block matrix with each block made up of matrices of size 2^^^. The initial matrix for 
L — 1 is given by 
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Equation is proved as follows. Let v,v' denote binary vectors of length L and w,w' denote vectors of length 
L — 1. The 2^ binary configurations of the system of size L are listed increasingly from 00 ... 00 to 11 ... 11. If 
and Vi are two binary configurations the (i, j)th entry of the Markov matrix Ml is the rate of the process Vj Vi for 
i ^ j and the diagonal entries are given by — J2j^ii^L)j,i- Decompose the Markov matrix of size 2^ in four blocks 
of size 2^~^ according to the first bit, 



, Moo Moi 
= ' Mio Mn 



(17) 



First let us consider Miq. This encodes all transitions v = Ow — > w' = Iw' . The only allowed transition is Ow ~^ Iw 
with rate a. Hence A/io — alL-i- Similarly, Mqo encodes transition of the form v = Ow v' = Ow' . If w ^ w' and 
the first bits of both w and w' are the same, the transition is either in the bulk or on the right boundary, which is 
encoded completely in the matrix for w — > m;', namely M^-i. There is, however, an additional transition which is 
present in the L — 1 system but is not present in the L system. Suppose that the first bit wi of w and the first bit w'l 
of w' satisfy wi = 1 — w[ and that all the remaining bits are the same, then w ^ w' with rate a because of ([2]). We 
therefore have to subtract these transitions from Afoo- Thus Mqq — M^-i — a{a ^ 1l-2)- Note that the sum of the 
elements of any of the first 2^^^ columns of AIl is indeed equal to zero, because in each column we have once added 
and once subtracted a and therefore the diagonal terms are unchanged. 

For Mqi, we consider transitions of the form v — Iw — s- u' = Ow' . If w = w' , this transition occurs with rate a. One 
other transition depends on the first bits of w, w' and assumes that all the other bits are the same. If wi = 1 — w'l, 
then w ^ w' with rate 1 in the system of size L because of the transitions ([T]). Thus Mqi = alL-i + (cr ® 1l-2)- 
Lastly Mil encodes transitions v = Iw ^ v' = Iw'. The argument is now similar to that of the Mqo case. If the first 
bits of w and w' are the same, then all the other transitions are encoded by Mi„i. Additionally, if wi — 1 — w'l and 
the other bits are the same, then there is a transition given by ^ for the system of size L—1 which is not present in 
the system of size L and we have to subtract this contribution. Besides, in Moi we have added 1 + a to each column 
and we have subtracted a in Mn, therefore, to ensure probability conservation, we must subtract 1 from the diagonal 
in Mil, which results in Mn = Ml_i — — Oi{(T ® 1l-2)- 

This system satisfies the transfer matrix Ansatz: one special solution can be constructed recursively as we now 
explain. In fact from the study of small systems (up to size L = 6), we found 2^ independent solutions of ([5]). This 
reflects the fact that the spectrum of M^ is fully embedded in that of Af^+i- 



Theorem 2 There exists a transfer matrix for the model which can be expressed by a recursion of order one. If one 
writes the transfer matrix from size L — 1 to size L by a block decomposition of matrices of size 2^^^ x 2^^^ as 



1,L 



(.-1) \ 
(L-1) . 



(18) 



then the matrix Tl^l^i can be written as 



T, 



L,L+1 



(L) 



(L) 



(19) 



with 



T.- 



(L-1) 



(L-1) 2r'^"^ 



This, along with the initial condition 



-Tr 

a 



(L-1) 



J 



and T2^'' = 



2T. 



(L-1) j,(L-l) 



'(cr® 1 



L-2J 



(L-1) 



(20) 



1+/3 + q;/3 a + /3 + a/3' 
a 1 



a + a/3 




a/3 
a 



(21) 



determines recursively a family of transfer matrices. 
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The proof is carried out by induction. We first suppose that equation ([5]) is satisfied for the transfer L — 1 ^ L. 
We then prove ^ for the transfer L ^ L + 1 using the conjectured formula ([^0]) for the transfer matrices. For 
convenience, we omit the subscripts denoting the dimension of the identity matrices. Using the decomposition (flSl) 
and IjlSp in the equation M^Tl^ix = Tl-i,lMl-i yields the following two identities: 



and 



Ml-iTi - a{a (g, l)Ti + aTa + (cr (E> 1)T2 = TiMl-i, 
aTi + ML-1T2 a{a (g) 1)T2 = T2ML-1. 



(22) 
(23) 



(Here, to simplify the notations we have written Ti instead of T^^~^\ T2 instead of and 1 instead of 1l_2)- 

We shall need two other identities satisfied by the transfer matrices defined by the recursion (HOI)- These identities 
will be proved by induction. 

The third identity is the following: 



aT2{a®l) = aTi - T2. 



(24) 



Assume that this is true for the blocks Ti,T2 of the transfer matrix Tl-i^ (fT5)) . Now, using (|^D|) . we have for the 
corresponding matrices in Tl,l+i 



which proves 

The fourth identity is 



(L) _ fTi+T2/a 2T2+T2ia(g,lL-2)/a 



T2 



T2I01 



2T2 T2(a®lL_2; 
T2 



T2(.^l,_2) 2T2\ ^^(.)(^, 



(25) 



^L-l 



T2[ML-l,{cTi 

This equation is proved by noting from (|15p that 

[Mi_i,(a® 1)] = 



1)] = [r2,(a® 1)] 



(cr® 1) 
-1 



1 

-(ct® 1; 



Then, using the induction hypothesis, we find that the left hand side of ((26|) is given by 



2T2 T2((7 1)\ fa(g>l 







T2 



-1 



1 



T2((7«)l) T2 

~T2 -T2{a®l) 



(26) 



(27) 



(28) 



which is easily verified to be equal to the right hand side of (j26p using again the induction hypothesis. 
Finally, using and (fT5|) . we calculate explicitly Tl^+iMl in square blocks of size 2^~^ as 



/ TiML-i+T2ML-i/a 



\ 

and similarly we write M^+iTl^l+i as 



aTi + 2T2ML-1 - 2T2 + Ti(cr ® 1) \ 



-ari(cr(g) l) + 2ar2 


-2aT2{a + T2{a (g> l)ML-i/a 


T2ML-1 + T2 


aT2 — T2/Q! 


-aT2(cr(g) 1) 


+T2ML-i/a 


2T2ML-i-aT2{(T®l) 


aT2 + T2{a ® 1) + T2{a ® 1)Ml-i 


aT2 


T2ML-1 - T2 ~ aT2{a (g> 1) ) 



(29) 



/ Ml-iTi + ML-iT2/a 
-a{a (g) l)Ti + 2Qr2 



aT2 



T2 + 2ML-1T2 - a{a (g> 1)T2 \ 
+ (cr (g> l)T2/a + Ml-iT2{(7 ® l)/a 
{a®l)T2{(T 



ML-IT2 


+ T2 


aT2 + ML-iT2/a 




—a{(T (g) 


1)T2 


-T2/a~ {a^l)T2+T2{cr<g 


)1) 


aTi + 2Ml- 


1T2 — T2 


2aT2 + ML-iT2{a®l) 




-2a{a(g 


1)T2 


+ (cr ® 1)T2 - a(cr (g) 1)T2{<7 d 


§ 1) 



(30) 



ML_iT2-T2-a(a(gl)T2 / 
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(Horizontal and vertical lines were inserted to separate the different blocks that compose the matrices). 

The fact that the matrices (j29l) and are equal is a direct consequence of the four identities (|22p , , (IM)) and 
(l26)l . This completes the proof of Theorem [2l 

There are many properties of the matrix T^^i+i which can be used to prove statements about the steady state 
distribution of the system. In particular, we can use (|12p and the recursion ([^0]) to compute the kernel The 

entries in vl are necessarily polynomials if we start with the vector \vi) — [ ^ | . Because the submatrix T2 of 

V " / . . 

is always proportional to a, it is easy to check that the value of for a configuration with k I's is always proportional 
to a'^ . In particular, the last entry of which corresponds to the configuration with all sites occupied is equal to . 

We shall call the sum Zl of the entries in vl, the partition function. This allows us to define the steady state 
probabilities as the vector \pl) = yz\'Vl)- Therefore the probability of the configuration with all sites occupied 

in a system of size L is given by /Z^. This corresponds to the smallest probability. Note that the polynomial 
Zl = ZL{a, (3) is the least common multiple of the denominators of the entries of the kernel pl of provided that 
the greatest common divisor of the numerators of the entries is one. For example, the system of size one has 

(31) 

whence Zi = a + 2(3. We find a remarkable property of the partition function of the system, namely its super-extensive 
growth with the size of the system. 

Corollary 3 The partition function of the system of size L is given by 

Zl = 2(^^'){1 + 2a)^-\l + p)^-\2a + 13). (32) 

This formula will be proved by induction. By definition of the partition function we have 

Zl+1 = (1l+i|ul+i) = {1l+i\Tl,l+i\vl), (33) 

where we have defined {1k\ = (1, • • • , 1| to be the line- vector of length 2^ with all entries equal to 1, K being an 
arbitrary integer. We now prove by induction that {Ik\ is a left-eigenvector of T^'^^ and T2^\ More precisely, we 
show that 




(U|t(^^= 2^-1(1 + a)(l+/3)(lK|, 
(U|Ti^)=2^-Vl+/3)(lK|. 



(34) 



For the initial condition in jUl), we have (1, l|rf ^ = (1 + «)(! + /9)(1, 1| and (1, 1\T^^^ = a(l + /3)(1, 1|. Then, using 
(pi)) we express T^'^'' in terms of T^'^^^K This leads to 

(Uiri^^ = 2(iA--i|ri^"'^ (35) 

(Here, we have used the fact that (lx-i| is an eigenvector of which in turn implies that {1k-i\T2^ = 

{lK-i\T^^~^\cr 1k-2))- Next, we use ([Ml) to obtain (l/flT^^^^ = (1 + '^/a){^K\T^^\ thus proving the recurrence 
([M]) . Finally, we have 

Zl+1 = (2^"'(1 + a)(l + (3){1l\ + 2^-1^(1 + (3){Il\)vl), 

= 2'^-\l + 2a){l+l3){lL\vL) (36) 
= 2^-1(1 + 2a)(H-/3)ZL, 



which, knowing Zi, proves the desired formula 

We now proceed by calculating some correlation functions in the model: 

Corollary 4 For a system of size greater than fc -I- 1, the probability of the first k sites being occupied is 

iVi--- Vk) = -7I^ ■ (37) 

2y^){l + 2a)'' 
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To prove this relation, we only have to consider a system of size + 1 as follows from Remark [TJ In this case, only 
two configurations contribute to this expectation value — either the last (i.e. the k + 1th) site is occupied or it is not. 
In the basis in which we have written our Markov matrix, this corresponds to the last two entries of \vk+i)- 

To obtain the sum of these two entries, we have to multiply the vector |wfc+i) on the left by the vector Wk+i = 
(0, . . . , 0, 1, 1| with 2'^+^ — 2 zeros. Using p^ . we have to calculate the action of {wk+i \ on the transfer matrix Tl-i^l. 
Note that w couples only with the last two rows of the transfer matrix. From the recursion (PO)) and the initial 
condition ((2T|) . we see that the only nonzero entries in the last two rows are from the 2x2 block at the bottom right. 
Notice also that the column- wise sum of these two blocks is the same and is equal to «(! + /3). Therefore, we obtain 



, , {wk+i\Tk+i.k\vk) a(l+/3)(0,...,0,l,l|i;fe) a{l + p)Zk{m ■ ■ -Vk-i) 
(r?i . . . r]k) = y = ^ ^ ^ 

^fe+l ^k+l 



(38) 



2'=-i(l + 2a) 

this recursion, along with the initial condition {rji) = yq^^, proves p7p . 

IV. CALCULATION OF ONE-POINT AND TWO-POINT CORRELATIONS 

A. Disorder Variables 

We now introduce a set of non-local variables that will allow us to calculate exact properties of the model. These 
are similar to the so-called interparticle distribution functions in [T^ . Consider the random variable 

(_l)#{i<j<»l';.=i} (39) 

which takes values ±1 depending on whether the number of occupied sites between the first site and the ith site is 
even or odd respectively. We call the disorder variable at site i. This definition is inspired from the theory of spin 
chains, where the Jordan- Wigner transformation allows us to write fermionic creation and annihilation operators from 
spin operators. More precisely, let Sj' , , Sj represent operators which raise, lower and measure the spin at site j 
respectively. They satisfy the relations {Sj'jS^} — 1 and [S'^,S'jr] — 0. From these operators one can construct 
fermionic operators fj , /J which satisfy the standard anticommutation relations 

/t = s;(-i)*>, 

/, = (-1)*<'>SJ, 

where = X]i=i(l/2 + S^) measures the number of up-spins to the left of the site j. Our variable £,i corresponds 
to the variable in the theory of spin chains. 

For an exclusion process like the one we consider here, has a more convenient representation, 

i 

c.=n(i-2^.)- (41) 

Note that these random variables satisfy £_f — 1 since the occupation variables satisfy 77^ = rji. The steady state 
expectation values of the variables satisfy some remarkable properties, which we will state in this section. These 
expectation values will be extremely useful in proving results for the density and the evaporation rate. 



Theorem 5 For a system with size L and i < L, 



an 



d 



{l3 + 2a){l + 2ay 



(a) = , , . • (43) 



9 



From (|¥T|) . we have for all i. 1 < i < L, 



from which we also have 



which substituting back in gives 



e; = (l-2r70e,-i, (44) 
Vi^i = -Viii-i, (45) 
+2?7,Cz- (46) 



We now write down the evolution equation for in the bulk. The main idea, and one of the reasons this variable 
is useful, is that a transition taking place strictly between 1 and i does not affect the value of (^i). The first bulk 
transition in ([1]) does not change the number of particles, and the second reduces it by two, and hence does not change 
the parity. We will first write down the equation and then explain each term. 

= -2a(6)-2(ry,6). 

The first two terms describe ways of entering the configuration contributing to If one starts with a configuration 
contributing to — (^i) and a particle enters or leaves site 1 the parity will change because of Similarly if one starts 
from configuration contributing to — (^i) and if the ith site is occupied and that particle leaves that site, then the 
parity changes. The last two terms describe ways of exiting the configuration contributing to (^i): when a particle 
enters from the left reservoir or when a particle occupying site i leaves it, we end up with a configuration contributing 
to — (Ci)- Using ([Tf|) in one ends up with 

(C.) = ^, (48) 
which along with the initial condition, (^o) = 1 leads us to (|42|) . For the last site, the balance equation is similar, 

- ^(a) = -2a(a) - WivUL) (49) 
because the last site exits with rate (3 instead of rate one. And using one obtains 

which leads to ([33]). 

We now calculate correlations among the disorder variables. These are directly related to the quantities of interest 
in the model. For example, 

= (1-27?,), (51) 

will give us the density. We will also see the relation between this problem and a random walk in two dimensions. 
For convenience, we introduce a new variable 

?m,n — 2 ' ^ ) 

and define the expression n+ to be the positive part of n. Namely n+ = n if n > and otherwise. 
Theorem 6 For a system with size L and < n < m < L, 

^ ' j=0 fe=l ^ ' 

and if L > n > m, then ^. 
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We first write down the balance equation for the correlation function for m,n> 1, 

(^^^"'^"^ 2(7ym^„j^,i) '^{'TlnS.m^n) J (54) 

where we have used ideas very similar to (j47p in the first line and (|46p in the second line. In this expression, is 
set identically to 1. The simplicity of the expression is due to the fact that the product ^mCn is affected only by 
transitions taking place between the mth and the 71th sites. Any other transition either changes the signs of both 
and ^„ (e.g, first site) or does not change the sign of either of them (e.g, last site). From equation ([51)) . we deduce 
the recursion relation 

;/■ \ _ (^m-l,rt) + {£,m,n-l) /^.p-x 
\t,m,n/ — ^ • (00) 



We now prove (|53p by showing that it verifies this recursion and the boundary conditions. Since the formula (j53p 
is valid in the triangular region < n < m, we have to check that it is valid for m = n and for n = 0. The case 
TO = n is easily done: substituting n = to in (j53p gives zero (because the sum is empty) and from (|52p . we also obtain 
S,m,m = 0. For the case n = 0, we have 

"2(^-2)+ (1 + 2a)" 2^1 h 

^ ' j=0 k=l ^ ^ 

/ \ (56) 

2'"-2l + 2a'»r ' ^\m-kl ' ^ 

^ ' \ fc=2 ^ ^ j=0 k=m-j 

where we have split the sum according to whether j = to — 1 or not, and changed the limits of k to count only the 
nonzero summands. The summation on k is easily done in both terms to give 

0) - « I + y\i + 2ay-A = (i+2«)--i 

I 1 + 2a ^ ^ ^ I 2(1 + 2a)™ ' ^ ^ 

which is the correct answer knowing ^^,0 = (1 ~ ?m)/2 and the expectation value (^m) from Theorem [SJ Thus we 
have verified the boundary conditions. We now need to verify the recursion (|54p . We can take to > n > 1 which 
implies to + n — 3 > 0. Then, 




(gm,n-l) + (Cm-1,«) _ « , r, xj-mr,(m-j-2)+ "^A^^ /to + n - 3 - (to - j - 2) + 

2 - 2.2™+"-3 + ^ ^ I TO-fc 

j=0 fc=l \ 

m — ri— 1 



— ^(l + 2a)J'-™+i2("-^-3)+ |^m + n-3-(TO-j-3)+ 



2 . 2m+«-3 A^"- I ^ \ m-k-1 

j=a k=i 



(58) 



We replace J J + 1 in the second sum to get the summands in both these terms to be exactly the same except for 
the binomial coefficients, in which the lower index differs by one. One can thus combine both these sums using the 
usual addition formula for binomial coefficients in the common range to get 

2 ^ 

j=l \ fe=l ^ ^ k=m-n ^ ^ / 

■m— n+1 

a 




, -(l + 2a)-™2'"-2 

2m+n-2 ^ ' / ' 

k=l 

(59) 

The second sum over fc, which comes from the contribution of (^m.n-i), involves only two terms which can again be 
summed using the addition formula for binomial coefficients and the result can be included in the first sum as the 
k = m — n term. The third sum over k (which is also an untouched term from (^m,n-i) corresponding to j = 0) is 
equal to one because the binomial coefficient is nonzero if n — 1 > to — fc, which is satisfied only for the upper limit, 
fc = TO — n + 1, where the value is one. Finally, this third sum over k adds to the first sum as the j = term, giving 
the right hand side of (j53p and confirming the recurrence relation. 



11 



B. Density 

We know from Remark [1] that the density at site k is a fixed quantity for sites L > k and depends only on a. From 
(|5ip . it is clear that rjm — S.m,m-i- Using Theorem [5] we arrive at an explicit expression for the expectation value of 
the density, 



- 2(2^-3)+ (1 + 2a)™ ^^^^ I m-1 



(60) 



When m = 1, we only have the j = term and the answer matches that given by Corollary 2] with k = 1. For m > 2, 
we split the sum according to whether j = m — 1 or not, 

We now multiply and divide this formula by two: the denominator in the prefactor becomes 2^™'~^ and the two terms 
inside the brackets are multiplied by two. For the first term, we get 

/2m -3\ /2m -3\ /2to - 3\ /2to - 2\ , , 

2(m-lj^(m-2j + (m-lj = (m-lj' ^''^ 

and for the second term, the power of two inside the summation becomes 2'""-'^^. We can then combine both these 
terms and write the density explicitly as 

= ——^ — y (^-^+ A (1 + 2a)^"2"-i-^". (63) 

22('»-i)(l + 2a)'" ^ V j J 

One can analyze the asymptotics of the density using Stirling's formula and noting that the largest summand is the 
one where j = m — 1. One finds that 

iVrn) - (64) 

Z^TTm 

independent of a. This asymptotic result was obtained by Lebowitz, Neuhauser and Ravishankar using stochastic 
couphng methods (see Theorem 2 of [IJ]). The problem studied by these authors corresponds to the case a = 1 and 
was inspired by earlier studies of the Toom model [32l | on the semi-infinite lattice. 

C. Evaporation Rate 

Since, unlike the TASEP, particles are not conserved in the bulk, a quantity like the current that is constant across 
all bonds does not exist in the present model. But we can define an evaporation rate at the bond (m — l,m), as 
the correlation {rjm-iVm}- this quantity does not depend on the size L of the system as long as i > m because of 
Remark [TJ We can calculate the evaporation rate using the expansion 

(Cm,m-2) = (1 - 277m " 277^-1 + "illm-lVm) , (65) 

from which we have 

(7?m-l77m) = • (66) 

An alternative method is to use the balance equation for the density 

^(?7m) = iVrn-l) - {Vm) ~ 2{rijn-lVm) = 0, (67) 

to obtain 

/ \ " ■r;^m-l-l/TO-2 + z\ ,v,„_2-i /^o\ 
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The asymptotics for the evaporation rate are given by 

iVm-iVm) - 7-7^= ■ (69) 
8v Trm'' 

One can similarly compute higher order correlations of a bunch of consecutive sites (77^ . . .77i+fe) using Theorem [S] 
The fact that the asymptotic evaporation rate in is given (upto a factor of 2) by the derivative of the asymptotic 
density (|64)l can probably be explained by looking at the hydrodynamic equation for this model. 

D. Structure of higher correlation functions 

We outline some general properties of correlation functions in the model. The reason one can make strong statements 
about correlations in this model is explained in Remark[TJ We will comment on correlations for both order and disorder 
variables. 

1. Density correlations 

We earlier used the balance equation for the density ([57)) in computing the evaporation rate. We also found simple 
closed form expressions for the correlations of sites 1, . . . , fc being occupied in Corollary [4| and for correlations of the 
form (11^=1(1 ~ '^Vj)) ill Theorem [5] In Theorem [S] we found a more complicated expression for correlations of the 

form (Jl^=/(1 ^ 277j)), where I 7^ 1 and I < k. Here, we will say more about the balance equations obeyed by more 
general order variables. 

Consider a general correlation function of order n, written as (yyp^ • • -rip^). In general we know no closed form 
expression for such a correlation function, but we can write down the balance equation satisfied by this object. We 
group the positions ii, - ■ - in according to blocks of consecutive sites. The total number of blocks is given by j and the 
length of the Ith block is denoted by A:/ + 1. The blocks are thus labelled as (77^^, • • • ,r]i-^^ki) up to (77^^, • • • ^rji-^^.). 
Note that neighbouring blocks are separated by at least one site and that the size of a block could be one. 

Theorem 7 The steady state equation satisfied by the correlation function of sites ii, • • • , ii + • • • ,ij, ■ ■ ■ ,ij + kj 
where ii > 1 is given by 

" (»?Pi • ■ ■ = XI ( n • ■ ■ 



;=i \ \m=l 



mi-i{^-2rii^Hi+i---Vii+ki i Y[ 77i„ . . . 77,^+fc„ 



(70) 



Vm=/+1 



The main idea is to consider what happens for a single block made of consecutive n sites i, . . . ,i + n — 1. The 
balance condition for such a block can be written easily, 

= ^(j?* • ■ -r/i+n-l) =('7i-l(l - ViH+l ■ ■■Vi+n-l) 

- iVi-iViVt+i ■ ■■Vi+n-i) - n{r^i ■ --Vi+n-i), (71) 
=(?7i_i(l - 2rii)rii+i . ..rji+n-i) 

- n{rii . ..Tji+n-i), 

because the configuration with all the sites i , . . . , i + n — 1 occupied can only be reached if the (i — l)th site is occupied, 
the ith site is empty and the remainder are occupied; we can exit the configuration if the (i — l)th site is occupied 
along with all the others from z to i + n — 1 with rate one, and in n different ways (each site between i and i + n — 1 
could jump) otherwise. The expression (|7ip is precisely ([70)1 for j = 1. Two blocks are separated by at least one 
site and this argument applies independently to each block. From each block I we get a contribution of /c; + 1 times 
{rjp-^ . . . rjp^ ) all of which add to give the factor n on the left hand side of ([70)) whereas on the right hand side the 
first factor rj^ of each block I is successively replaced by a factor ryij_i(l — 2?7iJ whereas the other blocks are left 
unchanged. 
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2. Disorder variables correlations 



A nontrivial observation about one-point and two-point correlation functions of the disorder variables is that their 
evolution equations (|i7)) and , supplemented by the recursion (gS]) are closed in the sense that they involve only 
other one-point and two-point disorder correlations respectively. In general, one would expect a hierarchy where the 
equations for lower correlations would necessarily involve higher correlations making the problem extremely difficult 
to solve for generic n-point correlations. As we will show below, this property of decoupling, which is very specific to 
this model, holds for all n. 

Theorem 8 In the steady state, the correlation function of n disorder variables {^p-^ . . . ) satisfies the equation 

n 

fe=l 

if n is even, and 

n 

{n + 2a){^p, ...^pj = X!(^Pi ■■■^Pk-iipk-iipk+i ---ipJ, (73) 
fc=i 

if n is odd. 

To prove these identities, one simply writes down the evolution equation, which for even n is given by 

d " 

^(Cpi ■ • ■ = ~2 ^ ] {£,pi ■ ■ ■ ^Pk-i^PkVpk^Pk+i ■ ■ ■ Cpn); (74) 

k=l 

and for odd n is given by 

d " 

'Tii^pi ■ ■ -^Pn) ~ ~2q!(^pj . . .^p„) — 2 ^ ] {£,pi ■ ■ ■ £,pk-i£,PkVpk^pk+i ■ ■ -Cpn) ■ (75) 

k=l 

One has to treat separately the even and the odd cases because a particle entering from the left reservoir makes a 
difference to the product S,p-i ■ ■ ■ Cp„ only if n is odd. The terms inside the sum itself are easily explained analogous to 
(|T7)l and ([M)) : each factor rjp^ occurs because the jump of the particle at site pk changes only the sign of £,p^ and keeps 
all the others intact. Then we use the recursion for each term in the sum to prove the result. We emphasize that 
no problem arises if any of the p^'s are consecutive. If Pk-i = Pk — ^ for a particular value of fc, a factor Cpfc_i — 1 
appears in the expectation value and drops out from the product. As an example, the correlation function for three 
consecutive sites satisfies: (3 -I- 2a)(^j^i+iCi+2) = (6-i6i+iCi+2> + (6) + (Ci+2)- 



V. SPECTRUM OF THE MARKOV MATRICES 



In this section, we mention some spectral properties of the Markov matrix of the asymmetric annihilation model. 
We observed that its characteristic polynomials factorize into linear factors and can be written explicitly. We first 
define the polynomials Al{x) and Bl{x) as 

\L/2] 

^Lix) = Yl (a; + 2fc)(V.'), 

LL/2J ^ ' 

Bl{x) = (.T + 2fc + 

k=0 

Conjecture 1 The characteristic polynomial Pl{x) of Ml is given by 

Pl{x) = Al{x)Al{x + 2a + I3)Bl{x + I3)Bl{x + 2a), (77) 
and successive ratios of characteristic polynomials are given by 

= Bl[x + 1)Bl{x + 2a + /3 + 1)Al[x + (3 + 1)Al{x + 2a + 1). (78) 

Pl{x) 
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If this is true, the Markov matrix has only 2L distinct eigenvalues. In particular, the negative of all the factors 
in the denominator Zl are roots of the characteristic polynomial; for example, the denominator Z3 according to 
Corollary [3] is 2(1 + 2a)^{l + Pf{2a + (3) whereas Psix) is given by 

x{x + 2){l3 + x + lf {2a + x + lf{2a + 2 + x + P){2a + x + (3). (79) 

We also conjecture that the Markov matrices are maximally undiagonalizable in the sense that each eigenvalue 
seems to have exactly one eigenvector independent of the number of times the eigenvalue appears (this implies that 
the Jordan blocks are of the maximum possible size). In other words, the degeneracy seems to be so strong that the 
Markov matrix Ml which is of size 2^ has only 2L eigenvectors. We can also write Pl{x) in a different manner. We 
define ni(m), for any integer m, as (—1) raised to the number of ones in the binary expansion of m. For example, 
7ii(l) = ni(2) = -l,7ii(3) = 1. 

Conjecture 2 Let us index the rows and columns from to 2^ — 1. The characteristic polynomial of Ml is then 
given by 

2^-1 

Pl{x) = n - i^^L)i^,^) " ani(z)). (80) 

i=0 

Equation ([77)) has been checked for systems of length < 7. The proof of this conjecture seems to be a nontrivial 
problem in determinant evaluation ^34] . 

We have been able to prove a much weaker result. Consider the matrix Ml — Ml — a{a 1) written in the block 
diagonal decomposition as 



ML-i-a{a(g)lL-2) 




{a(g)lL-2) 
Ml-1 - -a(cr(8) li 



(81) 



using Notice that the two diagonal blocks are Ml~i and Ml-i — Using one can see that Mi is 

upper triangular. Therefore Ml is also upper triangular for all L. Therefore the eigenvalues of Ml are simply the 
elements on the diagonal. Using the special nature of the diagonal blocks, one can easily prove by induction that the 
characteristic polynomial of Ml is given by 

P^{x) = Ylix + a + k)(''^') (x + a + (3 + k)i''^') . (82) 

k=0 

Note the close similarity with the characteristic polynomial of Ml in (|77|) . 

VI. DISCUSSION AND CONCLUSION 

In this work we have studied a nonequilibrium system on a finite size lattice with open boundaries in which 
particles diffuse and interact through hard-core exclusion and pairwise annihilation. The breaking of detailed balance 
in the bulk of the system is ensured by the asymmetric hopping rules and by the absence of pair-creation process. 
Besides, the difference in the chemical potentials of the left and the right reservoirs is also a source of nonequilibrium 
behaviour. The bulk dynamics is characterized by a single dimensionless parameter A which represents the ratio 
between evaporation and hopping rates ([1]). For A = 0, the model is identical to the totally asymmetric exclusion 
process with open boundaries Q. In the present work we have derived exact results for the case A = 1. The fact 
that A = 1 is a special point can be understood if one writes the Markov matrix Ml of this stochastic process as a 
non-Hermitian spin chain operator using the Pauli matrices [23 , [33 | . For general values of A we obtain 

L-l 

Ml = M,,»+i +R + L, (83) 



where 



M — Q+<^- -4-- ^qzqz 1 Ul^ _ 1 ^ ^ + ^ 

R - P(st-'-^] . (84) 
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We recall that the spin operators are given by: = (^"'^1,5' = (^^| and — ^ \ ■ For A = 1, we 

observe that all the interaction terms of the type SfSfj^^ disappear from the spin chain operator: this corresponds 
to the free fermion point. This absence of interaction gives an alternative explanation for the simplicity of the 
model at the special value A = 1. However, we remark that at this special point, the boundary terms give a 
nonlocal contribution because of lack of periodicity and therefore the usual strategies for diagonalizing the quadratic 
Hamiltonian on a periodic ring (i.e. Fourier Transform or Bogoliubov transformation) do not seem to apply here. The 
transfer matrix technique allows us to by-pass these difficulties. For A = the spin chain operator ([5^ represents the 
totally asymmetric exclusion process which is also exactly solvable. It would therefore be of interest to explore the 
integrability properties of this system for general values of A taking into account the boundary conditions. Another 
possible extension is to allow backward hopping of the particles (which corresponds to the partially asymmetric case) 
and to formulate the boundary conditions so that the model remains solvable. 

For this asymmetric annihilation process, we have been able to derive exact combinatorial expressions for correlation 
functions such as local densities and evaporation rates. The use of disorder variables has been helpful. However, the 
main tool that we have introduced is a recursion relation between systems of two consecutive sizes. This recursion is 
encoded in a semi-similarity operator (the transfer matrix Ansatz) that conjugates the Markov matrices for systems 
of sizes L and L + 1. The model studied here admits a transfer matrix Ansatz whereas the TASEP does not (as 
we explicitly checked on small systems). Conversely, the TASEP can be solved using a quadratic matrix product 
representation (which allows the calculation of physical observables such as the density and two-point correlation 
functions) whereas the steady state weights of the asymmetric annihilation process cannot be written easily as a 
simple matrix product (see Appendix ■ Indeed, the fact that the normalization Zl grows super-extensively as 
2L(L-i)/2 j^ggg l\'S2\i ) implies that the matrices for the stationary weights should involve L{L — l)/2 tensor products 
as in the case of multi-species exclusion processes [H, and the calculation of physical observables would be a true 
challenge. 

To summarize, the method for solving this model differs considerably from that used for the exclusion process. The 
transfer matrix Ansatz encodes in a particularly efficient way recurrence relations for the steady state probabilities 
and it allows one to deduce information about the steady state in a rather elementary manner. We believe that the 
existence of a transfer matrix Ansatz is rather general: it applies to the multi-species exclusion process on a ring [37l | 
and perhaps in a non-obvious manner to the ASEP with open boundaries. 

Of a more fundamental interest is the following question. The system is out of equilibrium and therefore its steady 
state violates detailed balance: this implies the existence of elementary currents between microscopic configurations. 
However, due to the evaporation of particles, there is no obvious way of defining a conserved current in the system. It 
would be interesting to find an observable that demonstrates at the macroscopic level the breaking of detailed balance. 

Finally, we remarked above that the normalization grows super-exponentially with the size of the system for 
A = 1 whereas for A = 0, logZ^ is extensive in L. Although Z^ has no direct physical interpretation, it grows with 
L roughly as the ratio of the most probable configuration to the least probable one. It would be of interest to see 
when this transition from exponential to super-exponential growth occurs as A varies and to interpret it physically. 
More generally, a challenging problem is to calculate for arbitrary values of A the scaling function /(A) defined by 
Zl ~ exp(/(A)L^). Related questions were addressed in rice-pile models by P. Pradhan and D. Dhar [38| and in the 
so-called Raise and Peel model by de Gier and collaborators |39l |. 
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APPENDIX A: RECURSIONS FOR THE STEADY STATE PROBABILITIES 



In this appendix, we give an algorithm that allows to calculate recursively the steady state weights of configurations 
in a system of size L knowing the ones for size L — I. These recursions are based on numerical observations and are 
conjectural. They are stated in this appendix in order to be compared with the recurrences for the TASEP. They 
could also be used as a tool to build a matrix product representation. 
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We define the pushing operator P that acts on a block of I's followed by a block of O's. P{VO ) is given by a linear 
combination of all possible ("'Jj*') binary words with j I's and k O's. The coefficient of a given word is a power of 2 
as follows: each move of the 1 which began at position j (i.e. the rightmost 1) contributes a factor of 2; each move 
of the 1 which began at position j — 1 contributes a factor of 1/2; each move of the 1 which began at position j — 2 
contributes a factor of 2 and so on. We give some examples, 

P(IOO) = (100) + 2(010) + 4(001), 
P(1100) = (1100) + 2(1010) + 4(1001) + (0110) + 2(0101) + (0011). ^ ' 

In the first example, there was only one 1 and we gained a factor of two for every push of the 1. In the second example, 
we gained a factor of two for each push of the rightmost 1 but lost a factor of two for each push of the other 1. 

We shall now construct an operator N which acts on words and generates a linear combination of words of size 
one less. The operator N will be defined in two stages: first we specify the action of N on special words that we 
call sub- configurations (which are building blocks for general words), then we explain how to extend this action on 
arbitrary words by concatenation. 

A sub-configuration is a word that contains a single block of empty sites. In a bulk sub- configuration, of the type 
10*^1 with A: > 0, the block is surrounded by an occupied site on the left and on the right. In a right (left:) sub- 
configuration, of the type 10*"' (0*^1) with fc > 0, there is a single occupied site on the left (right). The configuration 
with all empty sites will be treated separately. We shall distinguish three cases: 



• Bulk sub-configuration: 



Lfc/2J 

^^(10*^1) = J2 2^'^'(l^^'+^0'=-2jy (A2) 



Left sub-configuration: the relation involves a, 



L(fc-1)/2J [k/2] 

iV(0'=l)= J2 2^P(l2J+i0'=-i-2j)_^^ ^ 2^p(l2J0'^-2j)^ (A3) 

J=0 j=0 

• Right sub-configuration: the relation now involves /?, 

L(fc-1)/2J Lfc/2J 
7V(io'=)=/3 J2 2^P(l2J+i0'=~i-2^') + (l + /3) ^ 2'=-i-^P(l2J0'=-2^'). (A4) 

We are now in a position to construct the operator N for any configuration except the one with all empty sites. 
Any nonempty configuration can be written as r = o'^^l'^O'^i . . . I'^gfem^ jj-^ ^jjjg notation, a bulk sub-configuration is 
10*^*1 for i = 1, . . . ,TO — 1, the left sub-configuration is 0*^"1 and the right sub-configuration is 10*^™ assuming ko and 
km to be nonzero. If either of them is zero, the corresponding boundary sub-configuration does not exist. The idea 
is to use the algorithms for each of these sub-configurations as defined in (jA2[) . (|A3[) and (|A4[) and concatenate. The 
actual formula depends on whether fco is positive or not: 

m— 1 

{O'"'!) ® (l''-i © (lO'^'l)) ® I'^-i © (10'="), 

Nir)^{ (A5) 
a0 (l''-i © (10'=^^)) © l''"-i © (lO*^"), 

i=l 

if fco = 0, 

where the symbol © denotes concatenation and where we use the summation formulae for the sub-configurations 
defined previously. We emphasize that the size of the configuration is exactly one less than r. Note also that the 
right sub-configuration term A^(10'='") is empty if km = 0. For example, 

A^(OIO) = ((1) + a(0)) © (/3(1) + (1 + /3)(0)), 

= /3(ll) + a/3(01) + (l+/3)(10) + a(l + /3)(00). ^ ' 
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We have thus defined the operator N for all words except those consisting only of O's. The operator N encodes 
recursions which are precisely those given by the transfer matrix Ansatz as we checked on system of size < 7. 

The configuration with all empty sites requires a slightly different algorithm: this is the only case where the recursion 
is different from the one provided by the transfer matrix in (|20p . In fact, N{0^) involves all 2^"-*^ configurations of 
size L — I and the coefficients depend on whether these smaller configurations end in a one or a zero. We need here 
a related but different pushing operator, P', which wc define presently. P'{VO'^) is again a sum of all possible (■'^j'') 
configurations but the coefficient (which is again a power of two) is assigned differently. We first assign 2'"' to the 
configuration l-'O'"' and this time, we divide by two for every push of the rightmost 1, multiply by two for every push 
of the second-from-right 1, divide by two for the third- from-right 1 and so on. For example, 

P'(110) = 2(110) + (101) + 2(011) 
P'(1100) = 4(1100) + 2(1010) + (1001) + 4(0110) + 2(0101) + 4(0011). ^ ' 

(Note the similarities and differences between (IA1[) and (jA7p .) The recursion for the configuration with all zeros is 
given by 

/L(L~3)/2J L(i-2)/2j 

N{0^)=pi 2^Pil^3+^0^-^~^^)®l + a J2 2^^'(1^^0^-2-2j)^ ;L 

( L(i-3)/2j L(i-2)/2j ^^ (A8) 

+ (1 + /3) a Y 2^P'(l2j + l0^-3-2j) ^ Q ^ 2^~lp,^^2jQL^2-2j^ ^ Q 

We remark that the powers of two are governed by the two pushing operators P and P' . Besides, except for the 
exceptional term with all zeros, there is a prefactor f] if the configuration ends in one and a prefactor 1 + /3 if it ends 
in zero. 

Equations (jA5p and (jASp fully define the action of N on an arbitrary binary word. The operator N acting on a 
word of size L returns a linear combination of words of size L — 1. If C is a configuration of the system of size L, then 
its unnormalized weight W{C) is defined as 

W{C) = {C\vl) , (A9) 

where vl is the steady state vector constructed using the transfer matrix Ansatz (fTT|) . We claim that W{C) can be 
calculated knowing the weights of the configurations of size L — 1 from the following formula: 

W{C) = W{N{C)) . (AlO) 

Namely, the weight of C is the linear combination of the weights of the configurations of smaller size generated by 
applying the operator to C. For example, from (jA6p . we deduce 

VF(OIO) = PW{11) + aPW{01) + (1 + P)W{10) + a(l + P)W{00). (All) 

A challenging problem would be to find a matrix representation (or more generally an algebra) that embodies the 
reduction rules (jASp and (|A8|1 . We believe that such an algebra does exist as in the case of ASEP. However, we 
emphasize that we have bypassed the matrix product representation altogether thanks to the transfer matrix Ansatz 
defined in 
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